Azimuthal asymmetries in single polarized proton-proton Drell-Yan processes 



Zhun Lu/'2 Bo-Qiang Ma,3'^'B and Jiacai Zhu^ 

'Department of Physics, Southeast University, Nanjing 211189, China 
^ Departamento de Fisica, Universidad Tecnica Federico Santa Maria, and 
Centro Cientifico-Tecnologico de Valparaiso Casilla 110-V, Valparaiso, Chile 
School of Physics and State Key Laboratory of Nuclear Physics and Technology, 
Peking University, Beijing 100871, China 
Center for High Energy Physics, Peking University, Beijing 100871, China 

We study the azimuthal asymmetries in proton-proton Drell-Yan processes with one incident 
proton being transversely or longitudinally polarized. We consider particularly the asymmetries 
contributed by the leading-twist chiral-odd quark distributions. We analyze the asymmetries with 
sm{2(f> + (ps) and sm{2(f) — 4>s) modulations in transverse single polarized p^p Drell-Yan and sin 20 
asymmetries in longitudinal single polarized Drell-Yan at the Relativistic Heavy Ion Collider, 
the Japan Proton Accelerator Research Complex, E906 (Fermi Lab), and the Nuclotron-based Ion 
Collider Facility (Joint Institute for Nuclear Research). We show that the measurements of the 
asymmetries in those facilities can provide valuable information of the chiral-odd structure of the 
nucleon both in the valence and sea regions. 
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I. INTRODUCTION 

The single spin asymmetry (SSA) appearing in vari- 
ous high-energy scattering processes is among the 
most challenging issues of QCD spin physics. Large SSAs 
were observed experimentally in the process pp^ tt X 
@ two decades ago. Standard perturbative QCD based 
on collinear factorization to leading power oil/Q cannot 
explain these asymmetries Q. Many theoretical stud- 
ies 0-[l| have been proposed to explain the origin of such 
asymmetries. One standard approach is to assume the 
existence of parton distribution and/or fragmentation de- 
pending on intrinsic transverse momentum, by going be- 
yond the collinear picture. In this transverse momentum 
dependent (TMD) framework, novel structures of the nu- 
cleon emerge. For instance, due to the correlation of nu- 
cleon transverse spin S and quark transverse momentum 
kT, there can be an asymmetric distribution of unpolar- 
ized quarks in a transversely polarized proton [loj : 

/,/pt (x, kr) - fq/p-t (x, -kr) 



2 fix (x, krp)- 



M 



(1) 



Here or A^/^/pt is referred to as the Sivers func- 
tion [^,1^ , and has been applied to explain the SSAs ob- 
served in the process pp^ — > tt X . For a while the Sivers 
function was thought to be forbidden by the time-reversal 
invariance property of QCD [Hi] . However, model calcu- 
lations [l^l by Brodsky, Hwang and Schmidt show that 
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the Sivers effect can be allowed in the semi-inclusive 
deeply inelastic scattering (SIDIS) and Dell-Yan process 
at leading- twist level, due to the final/initial state in- 
teraction (FSI/ISI) between the struck quark and the 
target remnant. It was then realized that FSI/ISI can 
be accumulated into the Wilson lines (gauge-links) that 
are the key ingredients for a full gauge-invariant defini- 
tion fll, [ij of TMD distribution functions. This also 
leads to the prediction on the sign reversal of the Sivers 
functions in SIDIS and Drell-Yan [l^. For hadron pro- 
ductions in hadron- hadron collision (i.e.. Ha + Hb 
hi + h2 -\- X), the situation is more involved, as there 
are colored objects in both the initial state and the fi- 
nal state. The multiple FSI/ISI will generate process- 
dependent TMD distributions [T6l - [l9| which are different 
from those in SIDIS or Drell-Yan process. This is also 
viewed as the breakdown of the generalized TMD factor- 
ization in inclusive hadro-production of hadrons [l^l . 

Allowing naive- T-odd parton distributions encourages 
a lot of theoretical and experimental studies. Substantial 
SSAs contributed by the Sivers effect in SIDIS processes 
[2l| - |28j . with one colliding nucleon transversely polar- 
ized, have been measured by several experiments dur- 
ing recent years. The asymmetries are identified by the 
angular dependence s\n.(4>h — 't's), where (ph and (ps de- 
note, respectively, the azimuthal angles of the produced 
hadron and of the nucleon spin polarization, with re- 
spect to the lepton scattering plane. The data on the 
Sivers SSAs have been utilized by different groups p9l - 
|33| to extract the Sivers functions of the proton, on 
the basis of the TMD factorization js^, HH] . Those sets 
of parametrizations of the Sivers functions were applied 
to predict the Sivers SSA in various processes, such as 
the SIDIS at Jefferson Lab (JLab), and the Drell-Yan 
processes at the COmmon Muon Proton Apparatus for 
Structure and Spectroscopy, the Relativistic Heavy Ion 
Collider (RHIC), and the Polarized Antiproton eXper- 
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iment (PAX). Many planned measurements of SSAs in 
single polarized Drell-Yan processes at the established 
or planned hadron accelerators/colliders have been pro- 
posed. One of the main goals of these experiments is 
to test the sign change of the Sivers functions in SIDIS 
and Drell-Yan process [sO, llBl , as a crucial prediction of 
QCD dynamics. It is also worthwhile to mention that a 
sign mismatch for the fc^-moments of Sivers functions has 
been found when the authors of Ref. [s^l compared the 
functions extracted from SIDIS data and those extracted 
from p^p — ^ -kX data. 

The planned polarized Drell-Yan processes at (future) 
available facilities also provide great opportunities to 
investigate various spin and transverse momentum de- 
pendent (TMD) distributions. Besides the Sivers ef- 
fect, there are some other effects that may contribute to 
the azimuthal spin asymmetries at leading twist thereby 
could be measured in single polarized Drell-Yan pro- 
cesses. It is interesting to point out that all these leading- 
twist effects (except the Sivers effect) involve the chiral- 
odd parton distribution functions. For example, the fol- 
lowing combinations 

hi(E)h^, h^j.(E)h^, h^^^h^, (2) 

will lead to SSAs with sin(2(/) — (ps), sin(20 -I- (ps) and 
sin 2(/) angular dependences, respectively. Here and (f>s 
are the azimuthal angles of the dilepton pair and proton 
transverse spin with respect to the hadron plane, and 
we use the convention for the angle definition introduced 
in Ref. [38*1 . These types of the asymmetries arise from 
the coupling of two different chiral-odd parton distribu- 
tions. The coupling hi (g) hj^ was first introduced and 
analyzed in Ref. [39| as an alternative mechanism for 
SSA and a method of accessing the transversity distri- 
bution functions hi [13, |4l| . The key ingredient for these 
SSAs is the Boer-Mulders function hi [42i| , which is also a 
naive- T-odd TMD distribution function and provides the 
necessary phase required for SSA. In this paper, we will 



present a phenomenological analysis of these SSAs in the 
proton-proton Drell-Yan process contributed by various 
leading-twist chiral-odd distribution functions. We con- 
sider proton-proton induced polarized Drell-Yan process, 
since there are several hadron accelerators/colliders, such 
as RHIC, the Japan Proton Accelerator Research Com- 
plex (J- PARC), E906 at Fermi Lab, and the Nuclotron- 
based Ion Collider fAcility (NICA) at the Joint Institute 
for Nuclear Research (JINR), that can perform these ex- 
periments. Therefore the asymmetries at different ener- 
gies and kinematical regions can be analyzed and com- 
pared, which is important for obtaining the information 
of various chiral-odd distributions functions from exper- 
iments. 

The remaining content of the paper is organized as 
follows. In Sec. II, we briefly review the systematics of 
leading-twist chiral-odd TMD quark distributions, then 
give the expressions of the corresponding azimuthal an- 
gle weighted asymmetries A^'^'^'^~'^^\ and 
^^uf'^ in the framework of TMD factorization. We con- 
sider both the single longitudinally and transversely po- 
larized Drell-Yan processes. In Sec. Ill, we present the 
phenomenological predictions for single transverse spin 
asymmetry in p^p Drell-Yan process, and single longitu- 
dinal spin asymmetry in p~^p Drell-Yan process at RHIC, 
J-PARC, E906 and NICA. We conclude our paper in Sec. 
IV. 



II. SYSTEMATICS OF LEADING-TWIST 
CHIRAL-ODD DISTRIBUTIONS AND THEIR 
ROLES IN SSAS 

At leading twist, according to the hermiticity prop- 
erties of the fields and parity invariance, one may de- 
compose the TMD quark- quark correlation matrix of the 
nucleon as follows |4a - |45| 



$(a:,fcT) = ^|/i?i+-/iT 



llT ■ 



M 



Sl /^1L 



Sl giL — 

kx ■ St 
M 



kx ■ St 
M 



'IT 



giTj 75 i+ 



2M 



2M 



(3) 



Here n+ = (0, 1, Ot) is a lightlike vector expressed in the 
light-cone coordinates, in which an arbitrary four-vector 
a is written as {a^,a^,aT}, with — (a*^ ± a^)/\/2 
and Ot ~ (a^, a^). The eight functions on the right-hand 
side of Eq. ([3]) not only depend on longitudinal momen- 
tum fraction a;, but also on the intrinsic transverse mo- 
mentum of the quark kx- Therefore they are named as 
transverse momentum dependent (TMD) distributions. 



or alternatively the three-dimensional parton distribu- 
tion functions (SdPDFs) in momentum space. As the 
extensions of the usual Feynman distribution functions, 
BdPDFs enter the description of various semi-inclusive 
reactions and encode a wealth of new information on the 
nucleon structures that cannot be described merely by 
the leading-twist collinear picture. 
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FIG. 1: Left panel: definition of azimutlial angles in Refs. 
[39| and [63 |: right panel: definition of azimutlial angles in 
Ref. [13 . Using the replacements (j) — > —0 and (^s — 0, 

the definition in the left panel is transformed to that in the 
right panel. 



Each of these eight SdPDFs represents a special parton 
structure of the nucleon. Five of them, the Sivers func- 
tion /.j^, the Boer-Mulders function hj-, the pretzelosity 
hj~rp, the transversal helicity gix, and the longitudinal 
transversity hj^j^, vanish upon integrating fcy) over 
fc-r- Particularly, two BdPDFs, the Sivers function and 
the Boer-Mulders function are naive-T-odd distributions 
and account for the SSAs in various processes. Among 
the eight BdPDFs, hir, hj^, h^rp^ and hj^i^ are chirally 
odd, that is, they describe densities of the probed quarks 
with helicity flipped. Except , other three chiral-odd 
distribution are T-even. The relation between Hit given 
in Eq. ^ and hi is 



2M 



{x, fcf.). 



(4) 



Since hi naturally appears in the expression of related az- 
imuthal asymmetries, our discussion on the transversity 
in the rest of our paper is based on hi rather than hix- 
The distributions hi and h^ describe the densities of 
transversely polarized quarks inside a transversely polar- 
ized proton and an unpolarized proton, respectively. The 
distributions h^j, and h^j^ arise from double spin correla- 
tions in the parton distribution functions (PDFs), repre- 
senting the densities of transversely polarized quarks in a 



transversely (but in a different direction) polarized pro- 
ton and a longitudinally polarized proton, respectively. 

Because of the chiral-odd nature of hi, h^, h^j, and 
hj^j^, in high-energy processes they have to combine to- 
gether with another chiral-odd object, i.e., with the 
Collins fragmentation function in SIDIS, or with another 
chiral-odd distribution function in Drell-Yan, to mani- 
fest their effects. This makes them rather difficult to 
be probed experimentally. As a result, they are less 
known than the chiral-even distribution functions. Any- 
way, there are some efforts to extract transversity from 
SIDIS data [H, [ij, and Boer-Mulders function from 
SIDIS and Drell-Yan data For and h^^, 

there are extensive model calculations j53 - l60{ and some 
proposals to measure them in SIDIS and pp Drell-Yan 
processes. 

All the leading-twist chiral-odd parton distributions 
can be probed in single polarized proton-proton Drell- 
Yan processes: 

//^(Pi) +p{P2) -^i*{q)+X^ £{l) + -f X. (5) 

Here we assume that one proton (with momentum Pi) 
is polarized, and f or — >■ denotes its transverse polariza- 
tion or longitudinal polarization. In leading order, the 
dilepton pair is produced from the annihilation of the 
quark and antiquark from each proton. We denote the 
momenta of the annihilating partons from polarized pro- 
ton and unpolarized proton as ki and fc2, respectively. 
Then we can define the kinematical variables as 



9 = 

Xi = 

y = 



l + l' = {q\qT,q''), Q' = 

kl_ Q2 

^p+' 



2Pi-q 
2 \X2 



Xl 



2P2 ■ q P, 



(6) 



In the Drell-Yan process, if the transverse momen- 
tum of the dilepton qt is measured, we can apply the 
TMD factorization 34 , [s^ |6l| which is valid in the re- 
gion q^ ^ to write down the differential cross section 
of processes at leading order as ^] 



da 



(1 



cos 



dxi dx2 d^qr d^l I 
+ \St 



' pi 



b L — - — sm 2(pl'j^ ^ 



LU 



sm 



sin(2</) + <As)F|!^''*+^^) 



-sin(20-</.5)P^';}^'^-^^^ 



(7) 



Here (j) and (ps are the azimuthal angles of l± and St lepton £ in the center-of-mass system of the lepton pair, 
with respect to the normalized vector h = Qt/Qt, re- In Eq. ([7]) we only give the terms appearing in ([5]), and 
spectively; and dV, = d cos 6d(j) is the solid angle of the other terms do not contribute in our analysis below. We 
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FIG. 2: The light-cone diquark model results of xhi{x) (left panels), xhj^,j^ (x) (central panels), and x/i^^; {x) (right panels 
of valence u and d quarks at = 1 GeV'^, respectively. 



note that in the literature there are different definitions 
of angles and (f>s, as shown in Fig. [T] In this work 
we adopt the definition in Ref. [ssj. Also, we apply the 



so-called CoUins-Soper frame |62l | , in which the structure 
functions are expressed as 



J 



^sin(2<^-0s) 
TU 

.sin(20+0s) 
TU 



T-isin 2tp 



C[.hf: 
C 



C 



= c 



2{h ■ kiT)[2{h ■ kiT){h ■ fczr) - feiT • fc2T] - fciVC^- ' k2T) 



2M% 



2{h ■ kiT)[h ■ k2T) - fciT • k2T 7± 



(8) 
(9) 

(10) 
(11) 



In above equations we have used the notation 
C[w{k^TMT)f9\ =^e2 j d^k^T <fk2 



C2T 



X (5^^^ [qt — kiT — k2T) w(fclT, k2T) 
+ f{x,,klr)9''{X2,klr)\ . 

(12) 

Thus all structure functions depend on xi, X2 and qT — 
\qT\- 

As shown in Eq. ([7]), the structure functions 



psin(2</)-0s) 
TU 



and F'r 



,sin(20-|-0s) 



TU 



contribute to the cross sec- 
tion in the case in which one proton is transversely po- 
larized (denoted by subscript T), and will give rise to 



sin(2(/) — (f)s) and sin(2(/) + (f)s) angular dependences, re- 
spectively. The structure function F^^'^'^ contribute to 
the cross section in the case in which one proton is lon- 
gitudinally polarized (denoted by the subscript L), and 
will give rise to a sin 2(/) angular dependence. Therefore 
one can define the following azimuthal asymmetries 



A. 



sin(20-0s) 
TU 



{Xi, X2, Qt) 



Fl 



sin(20-0s) 



TU 



<sin(20-f 0s) / \ 
Atu {X1,X2, qT) 



^UU 

sin(20+0s) 



A sin 2(f)/ \ 

^LU X2, qT) 



Fr 



TU 



F\ 



F^ 
^uu 

sin 2(/) 



LU 



^UU 



(13) 
(14) 
(15) 
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Our definitions for tlie azimuthal asymmetries are similar 
to the analyzing power given in [39j and are different from 
the transverse momentum weighted asymmetries defined 
in [63j . For experimental measurement of the asymme- 
tries given in Eqs.([T3|. p4| . and (IT5|) . the polar angle 9 
of the lepton I should be identified. As a compensation, 
larger asymmetries could be measured. 

One can also express the cross section of the Drell-Yan 
process, depending on y and as 

dCT _ 1 do ^^^^ 

dy dQ'^ d^QT dfl s dxi dx2 d^qr dfl 
At the region ^ Q^, the following relations hold 

Q 



Xi 



X2 



(17) 



Therefore we can define the y-dependent and Q^- 
dependent SSAs as 



J d^qr dQ^^F^u 
Jd^qr dQ^^Fij^' 
Jd^qT dyF^ij 



(18) 



(19) 



Jd^qrdyF^^' 

where we have used the short notes P — T or L, and a = 
sin(20 ± 05) for P = T and a = sin 20 for P ^ L. The 
integrations in Eqs. ([TB)) and are performed according 
to kinematical cuts or experimental acceptances. 



III. PHENOMENOLOGICAL ANALYSIS OF 
AZIMUTHAL ASYMMETRIES AT RHIC, 
J-PARC, E906 AND NICA 

In this section we investigate the prospects of experi- 
mental measurements on the azimuthal asymmetries de- 
fined in the last section at various facilities that can con- 
duct single polarized proton-proton Drcll-Yan processes. 
The proton-proton Drell-Yan process involves the annihi- 
lation of a quark from one proton and a antiquark from 
another proton. In order to calculate Fpjj, one needs 
to know the distributions hi, h^rp, hj^j^, and of both 
the valence and sea quarks. Although there are some 
extractions of transversity and Boer-Mulders functions 
from SIDIS and Drell-Yan data, most of the chiral-odd 
parton distributions are not measured and less known, 
especially those of sea quarks. In order to estimate the 
azimuthal asymmetries in pp Drell-Yan processes, we ap- 
ply the following ansatz: 

• For the Boer-Mulders functions h^'^{x,k'^), we 
adopt the result extracted from the unpolarized 
pd [iJl and pp [65*1 Drell-Yan data in Ref. [49|, as 
there is parameterization for both valence and sea 
quarks with the following form: 



(^x , kf-p ) 



•(i-xfmx) 



1 



exp 



fr2 



(20) 



where the subscript "6m" stands for the Boer- 
Mulders functions, and q = u, d, u, and d. We have 
ignored the contributions from other flavors, since 
they are assumed to be small. We note that the 
possible range of parameters Hq allowed by the pos- 
itivity bound for hj^ can be described by the coeffi- 
cient uj, namely, that the substitutions Hq — >■ uHq 
for q — u, d and Hq — > -^Hq for q = u, d will 
not change the calculated cos 20 asymmetry (con- 
tributed by /i^' X hf^) in the unpolarized pd and 
Drell-Yan data. The range of u) given in Ref. ^40] 
is 0.48 < u! < 2.1 and uj — 1 corresponds to the 
central values of Hq. However, for the azimuthal 
asymmetries given in Eqs. ^TE\\ and (dH), the varia- 
tion of u will lead to the change of the magnitudes 
of the asymmetries, and will be considered in our 
calculations. 

• For the T-even distributions hi, hj^j,, and hj^^ of 
valence quarks, there are considerable model cal- 
culations. We will deploy the calculation from the 
light-cone quark-diquark model. In this model, the 
Melosh-Wigner rotation ^6§\ , which plays an impor- 
tant role to understand the proton spin puzzle (67j 
due to the relativistic effect of quark transversal 
motions, has been taken into account. In prac- 
tice, the light-cone quark-diquark model has been 
applied to calculate the helicity distributions [68| . 
the transversity distributions [gl, [tO] and other 
SdPDFs m, m, [lH, and related azimuthal spin 
asymmetries in SIDIS processes [zl, . 

The light-cone model results for the distributions 
h'\^ (x , k% ) , hj j ,^ ( x, k"^), and hfl^^{x,k'^) are given 
as [5H,[^,|6alzi 



1 



r{x,k'r)=[frix,ki)--ftix,k 



--ftix,kl)Wlix,kl), 
/''(:E,4)=-i/f (x,4)W^^(a:,4), 



W'six,k^T) 



(21) 
(22) 



where j = hi, h^^, h^^, respectively, and the su- 
perscript "w" is corresponding to the valence distri- 
butions. Wg^y{x,k'^) are the rotation factors for 
the scalar or axial vector spectator-diquark cases. 
Their explicit form are 



W'^'{x,kl\ 



<-(x,4) 



{xMi 



■mqf 



with 



Md = 



{xMd + rriqY 








{xMd +mqY 


~\~ k^ 


2Mn{xMd + 


ruq) 
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H~ 


m'^ + k^ ^ mjj + 




X ' 1- 


X 



(23) 
(24) 
(25) 

(26) 
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FIG. 3: Azimuthal asymmetries A^^'^'^"'^^' (left panels), A^^'^'*^'*^'' (central panels), and A^^^*^ (right panels) at RHIC 
collider experiments as functions of the rapidity y and the dilepton mass Q, respectively. The dashed lines in the left panels 
correspond to the contributions from the valence transversity distributions h\ fitted by Anselmino et al. in Ref. [13]. The thick 
solid lines in the left, central, and right panels represent the contributions from the distributions hi, h^rp, and h^j^ of valence 
quarks alone in the light-cone quark-diquark model [ssl . Is^ . Ig^ . ItoI ]. The shaded regions give the ranges of the asymmetries 
by considering the additional contribution from the distributions hi, h^j., and h^^^ of sea quarks constrained by the positivity 
bounds given in Eqs. (|27l) . (|28|) and (|29p . The upper and lower limits of the bands correspond to the asymmetries by saturating 
the positivity bounds. 



An important feature manifested by these rota- 
tion factors is that they automatically satisfy the 
requirement of the positivity bounds j74l] for the 
PDFs. In the left, central, and right panels of 
Fig. [2] we plot the curves for xhiix), xhf^^\x), 
and xhf^^'' (x) of valence u and d quarks at 
GeV^, respectively. 
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As there is already extraction of transversity from 
the global analysis by combining the SIDIS and 
e'^e~ annihilation data, we also use the most re- 
cent parametrizations [431 ^'^^ calculate the 
asymmetries A^'^^'^, and compare the results with 
those predicted from our model calculation. 



• In order to consider the effects of the distributions 
hi, h^rp and /i^^ of sea quarks, we constrain them 
by the positivity bounds [tJ 



h\{x, k^) ^ fi{x, fcfn), 



h-^^(x, k^) ^ fi{x, kx), 



2M% 



(27) 
(28) 
(29) 



They give rise to additional contributions to the 
asymmetries through the coupling (g) h^"^ , and 
will give a range of the asymmetries by varying the 
distributions within the bounds. By saturating the 
positivity bounds, one can obtain the upper and 
lower limits of the asymmetries. 



For the unpolarized distributions fi{x, fc^), we use 
the MSTW2008 LO set parametrization (tH, and 
adopt a Gaussian form factor for the transverse 
momentum dependence which has been adopted in 
many phenomenological analyses [i^ \^ 



f!{x,kl) = f!{x) 



(30) 



with kl^ = 0.25 GeV^; the subscript "wn" stands 
for the unpolarized distributions. 

• In order to precisely predict the azimuthal asym- 
metries at different experiments using TMD fac- 
torization, it is essential for one to know the evolu- 
tion of SdPDFs. Unlike the PDFs in the collinear 
factorization approach, whose evolution has been 
well established by the DGLAP equation, the Q^- 
dependence of 3dPDF is not fully understood yet. 
In our practical calculations we assume that the 
scale dependences of 3dPDF and the spin averaged 
distribution function /i are the same. The same 
assumption has been applied in some extractions 
of 3dPDF m, 113-113 • To what extent that this 
approximation is valid still needs to be studied. As 
the asymmetries we calculate are ratios, we expect 
that our assumption on the scale dependence are 
reasonable. 

Now we have all the ingredients for estimating the az- 
imuthal asymmetries in single polarized proton-proton 
Drell-Yan processes. In the following, we apply the above 
ansatz to present our predictions and phenomenological 
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FIG. 4: Azimuthal asymmetries yl^'^'^'^"'*^) (left panels) , A^i^'^'^+'^s* (central panels), anci Af^/''' (right panels) at RHIC 
fixed-target experiments. 



analysis for forthcoming experiments at RHIC, J-PARC, 
E906, and NICA. 

• RHIC 

The original proposal of Drell-Yan experiment at 
RHIC employs two proton beams to collide at 
= 200 GeV or 500 GeV But recently there 
is also a new proposal to conduct a fixed-target 
experiment at -y/s = 22 GeV f7T\. We estimate 
the asymmetries for both the collider and fixed- 
target modes at RHIC. The longitudinal proton 
beam will be run in the coming years at RHIC, af- 
ter that, Drell-Yan program with tranverse spin will 
be conducted. We choose the following kinematics 
for collider experiment at RHIC-STAR (Solenoidal 
Tracker at RHIC): 

^/J = 200 GeV, 4 GeV < Q < 9 GeV, 
< < 1 GeV, - 1< y < 2. 

We constrain the kinematical cut at the low trans- 
verse momentum region such that ^ where 
TMD factorization dominates. 

In the left panels of Fig. |3] we show the estimated 
azimuthal asymmetries A^r^'^^'^~'^^'^ as functions of 
the rapidity y and the dilepton mass Q, respec- 
tively. The difference between two linestyles in the 
left panels is that for the dashed lines we use the 
transversity distribution for valence u and d quarks 
from the parameterization in Ref. (33| , while for the 
solid lines we adopt the results from the light-cone 
quark-diquark model [tO] for the valence distri- 
butions hi. The thick solid and dashed lines cor- 
respond to the contribution merely from the com- 
binations of the valence transversity distributions 
and the sea Boer-Mulders distributions, that is, ig- 
noring the transversity distributions of sea quarks. 
The shaded regions give the ranges 



by considering the additional contribution of the 
transversity distributions of sea quarks constrained 
by positivity bound (|27p . The upper and lower lim- 
its of the bands correspond to the asymmetries by 
saturating the positivity bound. The first, second, 
and third rows show the results for lu = 0.5, 1 
and 2, respectively, where uj is the parameter for 
Boer-Mulders functions, as explained previously. 
In the central and right panels of Fig. [3] we show the 
estimated azimuthal asymmetries A^^'^"^^"^^' and 

^Lfj^"^' respectively, in the same way as in the left 
panels. 

From Fig. [3] we observe that in the forward rapid- 
ity region, the asymmetry A^'^ ^'^ at RHIC col- 
lider experiments is positive, while the asymmetries 
j^smj2<p+<ps) ^"^^"^ tend to be negative. It is in- 
teresting that the magnitudes of the asymmetries 
increase as the rapidity increases. At large forward 
rapidity, the asymmetry A^Jj/^*^ is dominated 
by the combination of the transversity of valence 
quarks and Boer-Mulders function of sea quarks 
(showed by the thick solid lines). This is under- 
standable since large y corresponds to larger xi 
and smaller X2. Therefore the measurement of the 
asymmetries at large rapidity can provide the in- 
formation of T-even chiral-odd distributions in va- 
lence region. Of course, the statistics at large rapid- 
ity are much lower than at midrapidity, therefore 
a reliable measurement requires data with high in- 
tegrated luminosity. A common feature shared by 
all these asymmetries is that as lu increases, the 
asymmetry in the forward rapidity region tends to 
decrease. This arises from the fact that larger uj 
corresponds to larger valence Boer-Mulders func- 
tion and smaller sea Boer-Mulders function. 

In the left panels of Fig. [31 the magnitudes of the 
asymmetries A^'^^^"^ calculated by using two 
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FIG. 5: Azimuthal asymmetries A^i^^^'^"'*^' (left panels) , ^5^^'^'^+'*^^ (central panels), and Afl,'^''' (right panels) at J-PARC. 



different forms of the valence transversity distri- 
butions hi are quite different at large rapidity re- 
gion. This is due to the fact that at large x region 
the valence transversity distributions hi fitted by 
Anselmino et al. in Ref. [l^ are smaller than the 
correspon ding ones in the light-cone quark-diquark 
model [iilTq]- 

There is also the possibility of accelerating the po- 
larized proton beam with Ep = 250 GeV to collide 
on the proton target at RHIC. The RHIC kinemat- 
ics for the fixed-target experiment are 

= 22 GeV, 4.5 GeV < Q < 8 GeV, 
< < 1 GeV, 0.2 < xi < 0.6, 

corresponding to —0.6 < y < 1.0, which is comple- 
mentary to the collider kinematics. In Fig. 2] we 
show the azimuthal asymmetries A^'^^'^ '^^^ (left 

panels), Af^.^j'^'^'^'^^^ (central panels), and A^l^^'^ 
(right panels) at RHIC fixed-target Drell-Yan pro- 
cesses. It seems that the magnitude of the asymme- 
tries j4^^^^"*~"*^^ and A^^'lf'^ at fixed-target experi- 
ments are larger than that at collider experiments. 
Therefore, there is a good chance to measure larger 
asymmetries at the fixed-target mode. The draw- 
back is that, at fixed-target experiments, the uncer- 
tainty from the T-even chiral-odd distributions hi, 
hirp, and h^j^ of sea quarks at the negative rapidity 
region is larger than that at collider experiments. 
In both modes the asymmetries are consistent to 
zero at the large backward region and their size 
increases with the increase of the rapidity. 

• J-PARC 

J-PARC might measure azimuthal asymmetries 
given in Eqs. (jl8ll9l) in single polarized Drell-Yan 
processes at Ep = 50 GeV [T^j, corresponding to 
y/s ~ IQ GeV. The kinematical cuts at J-PARC 



are 

4 GeV < Q < 5 GeV, < < 1 GeV, 
0.5 <xi< 0.9, 

corresponding to < y < 0.69. The estimated 
asymmetries Af^^^-^'\ Afj}''*+^'\ and A%^^ 
are shown in the left, central, and right panels of 
Fig. [5l The figure manifests that the asymmetry 
^sin(20-0s) positive, while the asymmetry J^^"^ 
is negative in all the allowed rapidity region. This 
feature can be seen even by considering the uncer- 
tainty at the midrapidity region from the distribu- 
tions hi and h^j^ of sea quarks. Larger asymmetries 
^si„(20-0s) ^sin20 predicted at J-PARC 

than those at RHIC. The asymmetry 

is much smaller than other two asymmetries, as at 

RHIC. 

• E906 

There is a new proposal to use proton beams from 
the main injector at Ep — 120 GeV to collide on the 
polarized proton target (NH3) by E906 Collabora- 
tion [7^ at Fermi Lab. The polarized dimuon Drell- 
Yan program at E906 might be applied to measure 
the asymmetries defined in Eqs. ((T5)) . and (|15p . 
The E906 kinematics are given as 

\/i=15GeV, 0.3<a;i<0.7, O.K 2:2 < 0.3, 
< (7T < 1 GeV, 4 GeV < Q < 7 GeV, 

corresponding to < y < 0.76. The calculated 
asymmetries Ap^'"^-'^'\ and A%^^ 

are shown in the left, central, and right panels of 
Fig.li 

• NICA 

NICA at JINR might realize both longitudinally 
and transversally polarized beams of protons at 
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FIG. 6: Azimuthal asymmetries A'^'^^^'''"'''^^ (left panels) , (central panels), and Afl^^''' (right panels) at E906. 



Vs = 12 27 GeV [80|. The kinematics applied in 
our calculation at NICA are 

yG^27 GeV, 4 GeV <Q <9 GeV, 
< gT < 1 GeV, 0.1 < xi < 0.8, 

corresponding to —1.1 < y < 1.1. Here we choose 
the highest cm. energy to avoid the overlap with 
other experiments. We present the asymmetries 



.sin(2<^- 



0s) 



A. 



i(2</,+0s) 



TU 



and A'21j^'^ in the left, cen- 



tral, and left panels of Fig. [3 respectively. 
Our predictions at RHIC, J-PARC, E906, and NICA 



show that the asymmetries A^^'^'^ "^^^ 



A 



sin(20+0s) 
TU 



, and 



^Lf/^*^ are sensitive to the Boer-Mulders functions of 
sea quarks. This can be seen by comparing the plots 
in the three rows of each figure. The size of sea con- 
tent can be described by the parameter uj appearing 
in the parametrizations of the Boer-Mulders functions. 
The case uj = 1 corresponds to the central values of the 
Boer- Mulders functions (which we refer as normal case), 
while w = 0.5 corresponds to much smaller valence and 
much larger sea values (large sea case), and ui = 2 cor- 
responds to much larger valence and much smaller sea 
values (small sea case) compared to the central values. 
In the normal case, we can see from the Q-dependent 
plots that the asymmetries A^^^^*^ and Af'^,'^^ are 



sizable at all entire allowed Q regions. As uj increases 
or decreases, the asymmetries decrease or increase corre- 
spondingly. Therefore, their measurements could be used 
to discriminate different scenarios of the Boer-Mulders 
functions. 

At larger backward rapidity region (RHIC and NICA) 
or midrapidity region (J-PARC and E906), the figures 
show that there are uncertainties contributed by the un- 
known sea content of hi, h^rp, and h^j^ allowed by the 
positivity bounds, especially in the small sea case. In 
some cases the uncertainties are so large that the sizes 
and signs of the asymmetries can not been determined. 



Precision measurement at these regions will provide fur- 
ther constraints on the sea content of hi, h^j,, and h^j^. 



The plots for Apf^-^'^ 



and A''11j^'^ show that these 



asymmetries are larger at the forward rapidity region, 
about 10% in magnitude in the normal case. Further- 
more, our plots show that at the forward rapidity region, 
the contributions from the T-even chiral-odd distribu- 
tions of valence quarks dominate, that is, they are less 
contaminated by their sea content. Hence the asymme- 
tries at forward rapidity are measurable and the measure- 
ments on them are ideal to access the valence content of 
hi, and /i^^ at large x region. 

The magnitudes of the asymmetries A^^*-^"^ "^^^ are 

larger than A^'fJ,^*^ and A^r^'^'^'^^'^'^K This is because that 
the size of the transversity distributions in the light-cone 
quark-diquark model is larger than that of h^rp and h^j^ . 
The comparison of different types of asymmetries might 
be used to distinguish the sizes of different T-even chiral- 
odd distributions and to check the approximate relations 
among TMDs [8l|. 



IV. CONCLUSION 

We have studied the azimuthal asymmetries in the sin- 
gle polarized proton-proton Drell-Yan processes by con- 
sidering particularly the contributions of the leading- 
twist chiral-odd quark distributions, i.e., the Boer- 
Mulders function, transversity, pretzlosity and longitu- 
dinal transversity. We define the azimuthal asymme- 

in transverse single po- 

sin 20 



tries ApJ^'^-^'^'^ and ApJ^* 



larized p^p Drell-Yan processes, and A^'j^"*" asymmetry 
in the longitudinal single polarized p^p Drell-Yan pro- 
cesses. Using the predictions for the transversity, pret- 
zlosity and longitudinal transversity from the light-cone 
quark-diquark model, and the Boer-Mulders functions 
extracted from the unpolarized Drell-Yan data at low 
transverse momentum, we present a comprehensive phe- 
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FIG. 7: Azimuthal asymmetries A^^'"'*'"^^^ (left panels) , yl^^'^'^+'^s) (central panels), and (right panels) at NICA. 



nonienological analysis of the asymmetries j4^^^^^~'~'^^\ 
^sin(20-0s)^ and J^^^^'^ at RHIC, J-PARC, E906, and 
NICA. In all these facilities there arc polarized Drell-Yan 
programs in preparation or being planned, including col- 
lider experiments (RHIC and NICA) and fixed-target ex- 
periments (RHIC, J-PARC, and E906). Our study shows 
that the polarized Drell-Yan programs at various facili- 
ties can be used to explore the valence and sea content 
of the leading-twist chiral-odd distributions in wide kine- 
matical regions. 
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